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New, high statistics measurement of the ir^e^u (i^es) branching ratio 

Julia A. Thompson, D.E. Kraus, and A. Sher*, for the E865 collaboration 
Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA 15260, USA 



70,000 K+ events to measure the K^^ branching ratio 



E865 at the Brookhaven National Laboratory AGS collected 
relative to the — > tt'^tv'^, — > tt^/x and K'^ — > Tr^Tr^Tr" decays. The tt" was detected using the e^e~ pair 
from tt" e^e~'y decay and no photons were required. Using the Particle Data Group branching ratios [0 for the 
normalization decays we obtain _Bii(_ft'J^^^j ) = (5.13 ± 0.02stat ±0.09sys ± 0.04,io, m)%, where ^s'jjj(^) includes the effect of 
virtual and real photons. This result is ^ 2.3a higher than the current Particle Data Group value. Implications for the 
Vus element of the CKM matrix, and the matrix's unitarity are discussed. 
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Figure 1. The cathedral at Durham. 



This contribution follows the thesis of Alexander Sher 
[|2j and the E865 group preprint [|3| prepared after 
the conference. I will comment on some differences in 
presentation, and emphasize systematic checks. Dis- 
cussions during all the conference activities (as for ex- 
ample in Figure^) were valuable in clarifying the con- 
nection of our data to other experiments. 

The experimentally determined Cabibbo-Kobayashi- 
Maskawa (CKM) matrix describes mixing between the 
" intrinsic" and physically observed quarks and is uni- 
tary within the Standard Model. One particularly in- 
teresting unitarity condition involves the first row el- 
ements: 



ud I 



IK. 



\Vub? = 1-5 



(1) 



where (5^0 indicates a problem with unitarity, and a 
sum less than one might indicate an additional quark 



generation. Since LEP results seem to have ruled 
out more than three neutrinos, such a shortfall would 
strain the standard model. Therefore, Eq. ^has gen- 
erated substantial interest. The Vud element is ob- 
tained from nuclear and neutron decays. While the 
precise shortfall in Eqn. varies with different deter- 
minations of Vud, the qualitative result is the same. 
For example, using an alternate, precise, recent value 
from the nuclear superallowed Fermi beta decays leads 
to (5 = (3.2 ± 1.4) • 10^3 [ 5,^ y^^^ too sjnah to affect 
Eqn. ^ is determined from semileptonic decays of B 
mesons 

Although Vus can be determined either from hyperon 
or from K — > T:eu decays, the purely vector _Ke3 de- 
cays, with less intrusion of hadronic physics, provide a 
smaller theoretical uncertainty [m^j. Theoretical con- 
tributions Vus were reevaluated recently [ 171 1^ ITnj . 
but since uncertainties of and \Vus\^ are com- 

parable, a high statistics measurement of the A'^ 
branching ratio (B.R.) with good control of system- 
atic errors is welcome. 

The bare (without QED corrections) ATij decay rate [ 



El [7| ISl ^3 can be expressed as 



dV{K+)^C{t)\Vus 



(2) 



'Present address: SCIPP UC Santa Cruz, Santa Cruz, CA 
95064. 



where t — {Pk — Ptt)^, C(t) is a known kinematic 
function, and /+(0) is the vector form factor value 
at t = which has to be determined theoretically [ 
l^'T'. Two recent experiments[ ll2l lll-t' provide A+ (the 
form factor slope) measurements consistent with each 
other and with previous measurements. An omitted 
negligible term in Eq. [21 containing the form factor /_ 
is proportional to M^/M^. 

E865 [E] searched for the lepton flavor number vio- 
lating decay — > 7r+/x+e^. The detector (Figure^ 
resided in a 6 GeV/c positive beam. The first dipole 
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Figure 2. Plan view of the E865 detector with a simulated 
—* -K^e^u decay followed by tt" — > e+e~7. 

magnet separated decay products by charge. The sec- 
ond magnet together with four multiwire proportional 
chambers (MWPCs: P1-P4) formed the spectrometer. 
The particle identification system consisted of thresh- 
old multichannel Cerenkov counters (CI and C2, each 
separated into left and right parts, for four indepen- 
dent volumes) filled with gaseous methane (7* « 30, 
effe ~ 0.98[|2]), an electromagnetic calorimeter, and 
a muon range system. The D and A scintillator ho- 
doscopes gave left /right and crude vertical position. 
The muon system was not used for this analysis. 

The decay 7r° was detected through the e'^e^ pair 
from the tt*' — > e^e~^ decay mode, with the 7 de- 
tected in some cases. In order to eliminate the un- 
certainty of the 7r° e+e~7 B.R. (2.7%), and to 
reduce experimental systematic uncertainty we used 
the other three major decay modes with a tt^ in the 
final state {K+ tt+tt^{K^2), K+ tt° ^i+iy{K^3), 
— > 7r+7r'^7r*'(K^Tr3)) for the normalization sample 
("Kdal"). 

The data were collected during a one-week ded- 
icated run in 1998, with beam intensity reduced to 
K, 10'^ if+/2.8 second pulse, a factor of 10 less than 
the standard K'^ —>■ 7r"'"/x"'"e~ running intensity, and 
trigger logic prepared especially for this measurement. 

The normalization (Kdal) and A'jg data were collected 
by the "ELER" trigger, which identified e+/e~ pairs 
(one left, one right) and required time coincidence 
of the four Cerenkov counters with at least one D- 
counter scintillator slat on each (left and right) side of 
the detector. The prescaled Cerenkov efficiency trig- 
ger, requiring only three out of four Cerenkov counters 
(no D-counter requirement), was used to measure the 
Cerenkov and D counter efficiencies. About 50 mil- 
lion triggers were accumulated, with « 37 million in 
the ELER trigger. The ELER events were about 3/4 



accidentals, often with a muon from high momentum 
beam particle decays K fiv or tt — > pLV making part 
of the Cerenkov portion of the trigger. 

Off-line event reconstruction, using the spectrome- 
ter only, required a three charged track vertex in- 
side the decay volume. The Cerenkov and D counter 
efficiencies were obtained from the Cerenkov effi- 
ciency triggers. The redundancy of the MWPCs (4 
planes/chamber), and track reconstruction was used 
to extract MWPC efficiencies. The absence of the 
electromagnetic calorimeter from the trigger allowed 
its efficiency determination from the data. Each effi- 
ciency was measured over its appropriate phase space. 

The Monte Carlo simulation used GEANT [15,. All 
relevant kaon decay modes were simulated, measured 
efficiencies were applied, and accidental detector hits 
(from the data) were added to the simulated events. 
The PDG value A+ = 0.0278 ±0.0019 [C] was used for 
the if ^3 simulation. 

The radiative corrections are calculated both outside 
and inside the 3-body i^es Dalitz plot: 

1. Outside: the K'^^^ (inner bremsstrahlung and 
structure dependent) decays which have photons 
hard enough to move the events outside outside 
of the 3-body K'^^ Dalitz plot boundary. These 
were explicitly simulated [^2; ^i'^ 0.5% of 
the ifij process. This 0.5% is present in our 
data, increasing our observed events from the 
"bare" K'^^, but subtracted when the "bare" K^^ 
B.R. is calculated. 

2. Inside: The difference between the observed 
events and the events from the bare i^Tij pro- 
cess. Here we used the procedure of Ref. [0. 
The observed events include four terms: a) the 
bare process (|aoP); b) the K^^^ process (|a7p; 
c) the virtual photons \ay\'^ correcting the bare 
process (small compared to the other correc- 
tions), and d) the interference of the bare process 
with the virtual photon corrections (2Reaoat,). 
This difference is —1.3%, i.e., the virtual pho- 
ton corrections and their interference with the 
bare process decrease the number of physi- 
cally observable events. The overall difference 
between the "observed, acceptance corrected" 
events and the events from the "bare" process 

is . ^observed, acceptance corrected — -^bare(l- ~1~ 

0.005 - 0.013). 

An overall correction of 1.0232 (the short distance 
enhancement) is also applied [ 1181 |H1 d- Following 
Cirigliano, et al.[Cj, we do not apply this to our data 
in calculating the B.R., but apply it in the calculation 
of Vus from the B.R.. 
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The difference, inside the DaUtz plot, between our de- 
tector's acceptance for the "bare" Ke3 process and for 
the observed process is an increase of 0.5%. Initially 
we expected the acceptance to decrease, since radia- 
tive corrections shift data to lower electron energies, 
where we lose events. However, we also have decreased 
acceptance at high electron energies, where the elec- 
tron energy distribution itself peaks. 

For the tt" e'^e~j decay, radiative corrections have 
been taken into account according to Ref. [.ISt- 

Tabled shows how the data were reduced by event se- 
lection. Selection criteria, common to K^^ and Kdal 



Table 1. Effect of data selection cuts on sample size (1998 
data). There is some overlap (which is accounted for in the 
analysis) in the final Ke3 and Kdal sample because of the 
third track selection criteria. The K — * tt+tt" acceptance 
is « 1.2%. Ke3 acceptance is « 0.7%, somewhat lower 
because of the lower average momentum in the Ke3 
decay. The acceptance can be approximately understood 
by taking a factor of three loss for each charged particle, 
30% for the C ambiguity cut, and approximately a factor 
of 2 for other quality cuts 



Selection Criterion 


Remaining 




Sample 




Size 


Online trigger (2 tracks) 


37,676K 


2 track vertex 


W9000K 


3 track vertex 


2852K 


CAmbiguity 


w 1900K 


Vertex, aperture cuts 


710 K 


Mee < 0.05 GeV/c^ 


644K 


Low mass e+e- particle ID 


626K 


Particle ID of third track (e/7r,//i) 




"Normaliz." (Kdal):7r7r,Kmu3,7r7r7r 


Ke3 


« 555K 


71K 



samples, included requirements for a good quality ver- 
tex, for the three tracks to cross the active parts of the 
detector, for the low (Mee < 0.05 GeV) mass e+e~ 
pair to be identified in the Cerenkov counters, and for 
the second positive track to have less than 3 GeV/c 
momentum. The momentum cut avoids any situation 
in which a mode other than Ke3 can have the sec- 
ond positive track from a normalizer mode properly 
satisfy the Ke3 criteria. A geometric Cerenkov ambi- 
guity cut rejected events where the Cerenkov counter 
response could not be unambiguously assigned to sep- 
arate tracks (w 30% loss of both Ke3 and Kdal). 

The K^^ sample was then selected by requiring the 
second positive track to be identified as e"*" in 2 of 
the 3 electron detectors: CI, C2, and the calorime- 
ter. Events entering the normalization sample had no 



response in at least one of the two Cerenkov coun- 
ters. These criteria minimized systematic uncertain- 
ties [0, but resulted in a small overlap of the ifij and 
Kdal samples, which was taken into account in the fi- 
nal result calculation. Final acceptances for the three 
charged particles differed by no more than 4% among 
the three normalization decays 

The final signal and normalization samples were 
71,204 and 558,186, respectively. Extensive compar- 
isons between quantities in the data and simulation 
were presented in the thesis [0 and the HEP preprint 

Contamination of the ifi^ sample by other de- 
cays occurred when tt^ or fi^ from normalization de- 
cays were misidentified as e+ from K^^, or as a result 
of tt" — > e+e~e+e~ decay. Care was taken with the 
PWC simulation so that track chisquares and vertex 
distributions agree between data and Monte Carlo for 
well-measured tracks. Decays (w 10% of tt's) are mod- 
elled in the simulation. Agreement of the standard 
vertex quality between data and simulation indicates 
successful modelling of the decays at the level required 
for our measurement. Systematic uncertainties were 
estimated by variation of the vertex quality cut. 

Total contamination of the Ke3 sample was estimated 
to be (2.49±0.05stat±0.32sys)%, with the systematic 
uncertainty caused by the simulation accuracy of the 
Cerenkov counters' response to 7r+ and /z"*". Contam- 
ination due to overlapping events was (0.25 ± 0.07)% 
and (0.12 ± 0.05)% of the selected normalization and 
signal samples respectively. Figure |31 shows energy de- 
posited in the calorimeter by the e+ from the selected 
-fCij sample. The contamination is manifest in the 
minimum ionization spike at 250 MeV. The small ex- 
cess of data in the spike agrees with our contamination 
uncertainty estimate. 

The final K^^ sample included «30% events with the 
fully reconstructed 7r°. We used these events as a con- 
sistency check but did not require photon detection in 
our main analysis. This eliminated an additional sys- 
tematic uncertainty from photon detection and recon- 
struction in the calorimeter. However, lack of a 7r° re- 
construction increased vulnerability to contamination 
from upstream decays and photon conversion. Up- 
stream decays whose photon produced pairs before the 
decay volume were suppressed by requiring the three 
track vertex to be more than two meters downstream 
of the decay volume entrance. In addition, the results 
obtained from the two independent samples (one with 
and one without the reconstructed) did not show 
a statistically significant discrepancy. The decay vol- 
ume was evacuated to about 10"^ nuclear interaction 
length, which suppressed beam pion interactions. 

The K^^ sample size gives a statistical precision of 
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Figure 3. Energy deposited in the calorimeter by the 
second positive track from the selected K^^ sample (e+ 
which is not from the low mass e'^e~ pair). No calorimeter 
information was used for the identification. Markers 
with errors represent data; the histogram is simulation. 

0.4%. The systematic error was determined from sta- 
bility of the result under variation of the reconstruc- 
tion procedure, selection criteria, detector efficiencies 
applied to the Monte Carlo, and subdivision (in pos- 
sibly biased distributions) of both signal and normal- 
ization samples . No significant correlations between 
different systematic uncertainties were observed. Sys- 
tematic uncertainties are summarized in Table |2 



in the momentum (Figure Q and spatial distributions 
[12]. The systematic error was determined by dividing 
ifij and Kdal events in roughly equal samples using 
the relevant parameter as a separator and observing 
the result variation [|2|. The sensitivity of the vertical 
spatial discrepancy to the MWPC alignment and of 
the momentum discrepancy to the spectrometer pa- 
rameters [12 indicate possible origins of these discrep- 
ancies. 




P(GeV/c) P(GeV/c) 

Figure 4. Reconstructed momentum of the e'^ from the 
low mass e+e~ pair for the selected K^^ and Kdal samples. 
Histograms represent Monte Carlo; points with errors are 
data. Monte Carlo to data ratios are shown on the right. 



Table 2. Systematic uncertainty sources and estimates 
of their respective contributions to the error of the final 
result. The total error was calculated as the sum of errors 
taken in quadrature. 



Source of systematic error 


Error esti- 




mate 


Magnetic field 


0.3% 


Vertex, quality 


0.6% 


Vertex position 


0.2% 


Cerenkov Ambig. 


0.3% 


Mee cut 


0.2% 


Aperture 


0.2% 


{n/fi)^ iden. 


0.04% 


MWPC efEc. 


0.2% 


D ctr. efRc. 


0.15% 


Cerenkov effic. 


0.3% 


Sample contam. 


0.3% 


Vertical distrib. 


0.8% 


1 er momentum distrib. 


1.3% 


K^^ trigger effic. 


0.1% 




0.1% 


Total uncer. 


1.8% 



The two largest contributions to the error come from 
the discrepancies [(2 between data and Monte Carlo 



In addition to the check using fully reconstructed 
7r°'s, as a second consistency check, we estimated the 
TT+TT+vr^ decay B.R. relative to the Kdal sam- 
ple. The result was (1.01 ± 0.02) x Rpdg, where 
RpDG is the prediction based on the PDG compila- 
tion, and the theoretical prediction [E| was used for 
the 7r° e+e~7 decay rate. The 2% error was de- 
termined by combining all relevant uncertainties in 
quadrature, and was dominated by the uncertainty 
in the prescale factor of the trigger used to collect 
iC+ TT+TT+TT" events. This result checks both the 
treatment of decays and the particle identification. 

A third check compared the K'^^ B.R. from 1998 and 
1997 data. The 1997 i^Tij data used a trigger requir- 
ing hits in the calorimeter, A and D-counters. That 
trigger neither allowed measurement of these detector 
efficiencies, nor of the K^.^ trigger efficiency. While 
we did not use the 1997 data for our final result, the 
1997 ir+3 B.R. was statistically consistent (within one 
sigma) with that from the 1998 data. This agreement 
is important since the momentum comparison in the 
1997 data looks qualitatively different from the 1998 
data[|2]. A preliminary reconstruction version was 
used for the 1997 data, without the final magnetic 
field tuning and detector realignment. Our intuition 
is that the discrepancies in decay product momenta 
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and spatial distributions, which dominate the system- 
atic uncertainties, reflect our imperfect knowledge of 
the magnetic field and detector positions but do not 
bias our result beyond our estimated systematic er- 
rors. Further refining of this knowledge was judged 
to be a Herculean task (the magnet, e.g., had been 
decommissioned) and was not undertaken. 

We estimated the form factor slope A+ from both 
1998 and 1997 K^^ data samples. We obtained: 
A+ = 0.0324 ± 0.0044,.tat for the 1998, and A+ = 
0.0290 ± 0.0044sfat for the 1997 data, both consistent 
with the current PDG fit. 

After subtraction of contamination in the Ke3 sample [ 
m, our result is BR{K+^^^)/{BR{K^2) + BR{K^3) + 
BR{K^3)) = 0.2002 ± 0.00085tat ± 0.0036^;,^, where 
ifi^^^^ includes all QED contributions (loops and inner 
bremsstrahlung). As noted above, for this result the 
tt" was detected using the e~^e^ pair from tt'^ e+e^7 
decay and no photons were required. 

Using the PDG fit [ values for the normaliza- 
tion branching ratios we infer BR{K^^^^^) = (5.13 ± 

0.02s(at ± 0.09sys ± 0.04„orm)% where the normaliza- 
tion error is determined by the PDG estimate of the 
normalization B.R. uncertainties. The PDG fit to the 
results to the previous decay experiments yields 
BR{K+ TT^e+v) = (4.87±0.06)% [[T], « 2.3f7 lower 
than our result. 

As discussed above, the total radiative correction was 
0.8%, yielding BR{K+ -> TT^e+Ubare) = (5.17 ± 
Qmstat ± 0.09,2,, ± 0.04„o™)%. This diff-ers slightly 
from the thesis (5.16%) due to recalculation with the 
2002 PDG, while the thesis used the 2001 PDG values. 

Using the current PDG value for Gp, the short- 
distance enhancement factor SEw[Mp, Mz) = 1.0232[ 
HE], and our result for the bare ifij decay rate we 
obtain |Ks/+(0)| = 0.2239 ± 0.0022rate ± 0.0007a+, 
which leads to |14,| = 0.2272±0.0023™te±0.0007A+± 
0.0018/^(0) if /+(0) = 0.9874±0.0084[|ni[7|. With this 
value of Vus and Vud from superallowed nuclear Fermi 
beta decays[|S], we obtain 5 = 0.0001 ± 0.0016. 

This result is consistent with CKM unitarity, but our 
ifij^ result is 5.3% higher than the PDG 2002 value, 
with a statistical error 0.4%, systematic error 1.8%, 
and normalization error 0.7%. 

We conclude: a) the Ke3 B.R. may be somewhat 
higher than that listed in the PDG tables; and b) Ke3 
experiments are now precise enough that radiative cor- 
rections should be applied consistently to all entries in 
the PDG average. 

Even without our result, the V„s extracted from 
the PDG average iC^ decay rate is higher than the 
Vus from the rate. [ UHl g]. Our K^^ resuh 



increases this difference. K^^ decay measurements 
(both charged and neutral) in progress (CMD2, NA48, 
KLOE)[ 4 should clarify the experimental situation. 
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